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We perform a theoretical and experimental study of a system of two ultracold atoms with tunable 
interaction in an elongated trapping potential. We show that the coupling of center-of-mass and 
relative motion due to an anharmonicity of the trapping potential leads to a coherent coupling of 
a state of an unbound atom pair and a molecule with a center of mass excitation. By performing 
the experiment with exactly two particles we exclude three-body losses and can therefore directly 
observe coherent molecule formation. We find quantitative agreement between our theory of inelastic 
confinement-induced resonances and the experimental results. This shows that the effects of center- 
of-mass to relative motion coupling can have a significant impact on the physics of quasi- ID quantum 
systems. 



A key question in condensed matter physics is how 
the dimensionality of a quantum system determines its 
physical properties. Especially in one dimension the in- 
creased role of quantum fluctuations leads to the appear- 
ance of interesting phenomena which cannot be observed 
in higher-dimensional systems. This poses the interest- 
ing question of how to experimentally realize such one- 
dimensional (ID) systems in a three-dimensional (3D) 
world. This can be achieved by confining particles in 
a strongly anisotropic potential whose lowest transver- 
sal excitation is much larger than all other relevant en- 
ergy scales of the system. In this case a 3D system can 
be mapped onto a true ID system obtaining an effec- 
tive ID coupling constant #id which depends on the 3D 
scattering length a [1 . In such anisotropic confinement, 
ultracold atoms have been used to study, e.g., the Tonks- 
Girardeau [2H3] and super- Tonks- Girardeau [5] gas as 
well as the fundamental question of what constitutes an 
integrable quantum system [6]. 

Such experiments |5j [TUTU] often rely on the fact that 
it is possible to control the effective ID coupling strength 
#id by tuning the scattering length a with a Feshbach res- 
onance [11]. For a specific ratio of the scattering length 
and the transversal confinement length gm diverges 
to ±oo at a confinement-induced resonance (CIR) pQ. 
To distinguish these resonances in the elastic scatter- 
ing channel from the molecule-formation resonances we 
study in this paper we will refer to them as elastic CIRs. 

A common experimental approach [12 to character- 
ize such resonances has been to look for increased loss 
of atoms caused by enhanced three-body recombination 
in the vicinity of the resonance. However, this interpre- 
tation of the observed losses has been called into ques- 
tion by a recent experiment which observed a splitting 
of loss features under transversally anisotropic confine- 
ment [12 , although later theoretical works showed that 



no such splitting of elastic CIR can occur p~3j [14] . One 
proposed explanation for the splitting is based on the fact 
that the trapping potentials used in experiments are not 
perfectly harmonic. This leads to a coupling of center-of- 
mass (COM) and relative (REL) motion [T5HT7] , which 
in turn can lead to a coupling of two atoms in the ground 
state of the trap to a weakly bound molecular state with a 
COM excitation [18] (further elaborated on in [19 ). The 
occupation of the bound state is only possible because 
the excess binding energy can be transferred into COM 
excitation energy due the anharmonictiy of the confin- 
ing potential. This redistribution of binding energy to 
kinetic energy is an inelastic process and thus we refer to 
these COM-REL coupling resonances as inelastic CIR. 

In a many-body system losses at the inelastic CIR can 
be described as a two-step process: First, two atoms 
coherently couple to the COM-excited molecular state. 
Then, this molecule collides either with another molecule 
or an unbound atom, which leads to a deexciation of the 
molecule into a deeply bound state and subsequent loss 
of the involved particles from the trap. 

However, also different theoretical models have been 
developed to explain the observed splitting of the loss 
features in [12]. These argue with enhanced three-body 
effects in the vicinity of elastic CIR. One is based on mul- 
tichannel effects [20] , others [12] [21] on a Feshbach-type 
mechanism. In a many-body system as used in [12] differ- 
ent loss mechanisms are in principle possible and cannot 
be clearly distinguished by the experiment. A straight- 
forward, yet experimentally challenging solution to this 
problem is to eliminate three-body effects by investigat- 
ing a pure 2-body system. In this work we provide a di- 
rect experimental confirmation of the theory developed in 
[T8] by performing a theoretical and experimental study 
of two 6 Li atoms in an elongated trapping potential with 
a slight ellipticity. Ab initio calculations of the coupling 
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strengths, the widths and the positions of the coherent 
molecule formation at the inelastic CIR are found to be 
in quantitative agreement with the experimental results. 

To prepare a quasi- ID two-body system we follow 
the same preparation scheme as described in [lOj [22] , 
which has a fidelity of about 90%. The two particles 
are trapped in the ground state of a cigar-shaped po- 
tential with a mean transversal confinement length of 
d± = 0.486 ±0.006 jam and an aspect ratio of about 10:1, 
which is well in the quasi-lD regime [23]. The shape and 
anharmonicity of the potential have been characterized 
by precise measurements of the transition frequencies for 
exciting a single particle into the first and second excited 
level in the longitudinal and both transversal directions 

This two-body system is in absolute coordinates de- 
scribed by the Hamiltonian 

tf(r 1; r 2 ) = Ti(n) + T 2 (r 2 ) + Vi(n) + V 2 (r 2 ) + UQn - 

(1) 

where T\, T2, Vi, V2 denote the kinetic energies and po- 
tential energies due to the trap of particles one and two, 
respectively, and U the interatomic interaction. It has 
been demonstrated that sextic potentials, i.e. expansions 
of a sin 2 optical-lattice potential up to order six, 

V (r)= E Y5 V ^ f -l V ^ j4 + V ^ j2 ' (2) 

j=x,y,z 

are well suited to describe anharmonicity induced COM- 
REL coupling in single- well potentials [25], like the one 
used in our experiment. 

The stationary Schrodinger equation for the Hamilto- 
nian ([T]) can be solved exactly by the computational ap- 
proach described in [26] . Herein, the interaction potential 
is treated by a numerically given Born-Oppenheimer po- 
tential curve of a 6 Li system. The variation of the scat- 
tering length due to the magnetic Feshbach resonance 
can be modeled computationally by modifying the inner 
wall of the potential curve which effectively changes the 
scattering length of the system to arbitrary values [27] . 

For a two-particle system it is convenient to transform 
the Hamiltonian in REL and COM coordinates, r = 1*1 — 
?2 and R = \{t-\_ + T2), respectively, 

i?(r,R) = T REL (r) +T CO m(R) +V REL (r) 

+ Fcom(R) + U int (r) + W(r, R). (3) 

Vrel(f) and Vcom(R) are the separable parts of the 
sextic potential [25]. Thus, W(r, R) contains only the 
non-separable terms, i.e. a polynomial in r 2 R 2 , r 2 R^ and 
r^R 2 . The potential parameters Vj and kj are obtained 
by fitting the eigenenergies of a single particle in a sex- 
tic potential to the experimentally measured transition 
energies of a single particle in the trap. The fit results 
are given in [24] . The eigenenergies and wavefunctions of 



the Hamiltonian (|3| can now be calculated via exact di- 
agonalization for different values of the s-wave scattering 
length. A fully coupled spectrum is shown in Fig. [I] 




Figure 1. (color online) Eigenenergy spectrum of the Hamilto- 
nian ([3]) for 6 Li atoms confined in a sextic trapping potential. 
In the upper part all states bending down to — 00 are molec- 
ular states originating from the REL bound state ip^ with 
different COM excitations. The two bound states marked in 
red are the only ones which have a significant coupling to the 
repulsive state (blue) . For the other states (gray) the coupling 
is negligible. The magnified part shows the avoided crossings 
responsible for the COM-REL resonances. 

Relative motion bound states ip^ with COM excita- 
tion <£ n , n = (n xi n yi n z ) (i.e. states bending down to 
negative infinity) cross with trap states, i.e. states whose 
energy converges asymptotically to a constant value for 
a — ^ + . In the absence of a trapping potential these 
states would lie in the continuum. The system is initially 
in the lowest trap state, i.e. dominantly in the relative 
motion repulsive state ^0 and COM ground state $ (0,0,0) 
(see blue state in Fig. [I]). Hence, it suffices to consider 
crossings with this state. The coupling, and equivalently 
the size of the avoided crossings, is described by the cou- 
pling matrix elements 

W n = <^ (b) $ n |W^o$(o,o,o)>- (4) 
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In [18] it was demonstrated that in quasi ID only 
the lowest transversally COM excited bound states, 
|^ (6) $(2,0,0)) and $0,2,0) (see red states in Fig. [I]), 
couple significantly with the lowest trap state. There- 
fore, in Fig.[l]only the transversally excited bound states 
form significant avoided crossings with the repulsive trap 
state. Due to the transverse anisotropy of the trap these 
crossings are non-degenerate which results in a splitting 
of the resonances. Such a splitting was also observed in 
[T2] in quantitative agreement with the positions of the 
inelastic CIR [18]. 

To demonstrate that the crossing states possess char- 
acteristics of a bound and a trap state the mean radial 
density 

POO 

f= drrp(r). (5) 
Jo 

was calculated. Here, 

p(r) = r 2 J dF R d^ r |^(r,R)| 2 (6) 

is the radial pair density where \I/(r,R) denotes the full 
six-dimensional wavefunction of the system, dVn is the 
COM volume element and dQ r is the angular volume ele- 
ment of the REL motion. At d y /a = 1.38 the bound state 
has a mean radial distance of f = 0.29 dj_, i.e. it is small 
compared to the mean transversal confinement length. 
This demonstrates the strong binding of the atoms. In 
the trap state, the atoms possess a mean distance of 
r = 1.19 = 3.06 d±. This mean distance which is of 
the order of the longitudinal trap length d z = 1.25 jam is 
a consequence of the elongated trap. 

In the vicinity of the avoided crossing the system can 
be approximately described as a two-level system because 
the other states are energetically almost inaccessible. 
When the scattering length is ramped non-adiabatically 
towards the crossing and stopped in the gap region of the 
avoided crossing the system finds itself in a coherent su- 
perposition of the two adiabatic states [28] , which in our 
case are the bound state <l> n ) and the repulsive trap 
state \ipo $(o,o,o))- Since both states evolve with different 
phase a Rabi-oscillation between the states occurs with 
the frequency 

n = i ^WITp (7) 

which is a measure for the coupling strength for S = 
(Eb—Et)/2 — 0. Here, W n is the coupling matrix element 
from Eq. Q while E^ and E t denote the energies of the 
diabatic bound and trap states, respectively. 

Experimentally this coherent superposition is realized 
by first preparing two 6 Li atoms in the ground state of 
the potential and then increasing the scattering length a 
by ramping up the magnetic offset field non-adiabatically 



with a speed of 20G/ms. To locate the molecule forma- 
tion resonances the ramp is suddenly stopped at differ- 
ent values of the magnetic offset field. The population 
is expected to oscillate between the unbound and the 
COM-excited molecular state as a function of the Rabi 
frequency Q which depends on the magnetic field. 

In a first experiment we wait for a fixed hold time 
of 12.5 ms after stopping the ramp at different magnetic 
field values between 779 G and 788 G [29]. We then mea- 
sure the number of free atoms remaining in the ground 
state of the trap by ramping to a magnetic field of 523 G 
where the molecules are deeply bound and therefore not 
observed with our detection scheme. Thus, the mean 
number of molecules is given by N mo \ = (Nq — Nqs)/Nq, 
where Nq is the mean number of atoms in the initial sys- 
tem and Nqs is the mean number of particles detected 
in the non-molecular ground state at the end of the ex- 
periment. To check whether the missing atoms indeed 
end up in the molecular state we repeated the experi- 
ment but ramped the magnetic field to a value of 900 G 
before measuring the number of particles. At this mag- 
netic field we are far above the elastic CIR so that the 
molecules become weakly bound and the constituent par- 
ticles of the molecules can be detected with our detection 
scheme. We found that there is no measurable change 
compared to the initial particle number when measuring 
above the elastic CIR, which excludes the presence of any 
significant loss channels in our system. Figure [2] shows 
the detected number of particles in the repulsive state de- 
pending on the magnetic offset field. As expected from 
numerics, two peaks are observable which are identified 
as the COM-REL motion coupling resonances created by 
the two molecular states excited in x- and y- direction of 
the anisotropic confinement. 

To analyze the dynamics of the coupling we ramped to 
different values of the magnetic offset field around the 
features shown in Fig. [2] and held the system for different 
hold times. With less than 10% probability we detect 
only a single atom in the trap, i.e. the two atoms are 
either free (N = 2 detected) or bound to a molecule 
(N = 0). The few realizations with just a single atom 
detected (N = 1) are not considered in the analysis. 
Fig.[3]a) shows the result of one of these measurements. 
The oscillation of the fraction of molecules shows that 
we have created a coherent superposition of the molecu- 
lar and the repulsive state. By performing a sinusoidal 
fit to the data we can extract the Rabi frequency ft of 
the oscillation. The maximum amplitudes of the oscilla- 
tion for different magnetic fields are shown in figure [3] b). 
From a Lorentzian fit to the amplitude we can extract the 
width (FWHM) of the coupling in terms of the magnetic 
offset field. Table [I] shows the width of the coupling res- 
onances determined from the measurement. The spacing 
between the two resonances, see Table [I] is in agreement 
up to 0.3 G with the numerical calculation. The absolute 
position of the experimental resonances is shifted about 
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Figure 2. Disappearance of particles in the repulsive non- 
bound state. Due to the COM-REL motion coupling the par- 
ticles in the non-bound state couple into a molecule and disap- 
pear when detecting the number of particles in the non-bound 
state. One observes two peaks indicating COM-REL motion 
coupling resonances involving two excited molecular states in 
x- and ^/-direction of the confinement. Each data point is the 
average of about 50 individual measurements with discrete 
atom number. The blue dashed line indicates the position of 
the elastic CIR at 779.3 ±0.5 G calculated using the transver- 
sal confinement length d± and the calibration of the scattering 
length a(B) of [30 as inputs for the theory of |31j . 



COM 


Position [G] 


FWHM[G] 


Q /2tt 


Hz] 


excitation 


exp. 


num. 


exp. 


num. 


exp. 


num. 


(0,2,0) 


780.5 


776.01 


0.25(0.03) 


0.35 


83(2) 


64 


(2,0,0) 


783.2 


779.02 


0.42(0.06) ( * } 


0.35 


75(1) ( * } 


69 



Table I. Comparison between experiment and numerical cal- 
culation. See [m for these measurements. 



4.3 G compared to the theoretical values. In view of the 
width of the elastic CIR [32 of 250 G this is a remarkable 
accuracy. Moreover, except for the two COM-REL cou- 
pling resonances no significant molecule formation was 
observed over the whole width of the elastic CIR. 

In conclusion, our results directly show that in a two- 
particle system the COM-REL coupling allows for the co- 
herent coupling of an unbound atomic pair and a molecu- 
lar state without a third particle being present. Compet- 
ing processes such as three-body recombination or pro- 
cesses involving atoms in higher bands [20 are excluded 
by our high preparation fidelity. Hence, the agreement 
between the theoretical and experimental results gives a 
quantitative confirmation of the theory of inelastic CIR 
[T8] . Furthermore, our results show that a molecule for- 
mation in a two-body system is absent at the elastic CIR 
[33] . The results strongly imply that the inelastic res- 
onances are the dominant cause for the appearance of 
the two distinct loss features in the experiment by E. 
Haller et al. [12] as was already suggested in [18]. In gen- 
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Figure 3. Coherent dynamic of the COM-REL motion cou- 
pling, (a) Oscillation between the non-bound and the COM 
excited molecular state. From a sinusoidal fit we deduce the 
Rabi- frequency Q. (b) Maximum amplitude of the oscillation. 
The data points are extracted from measurements analog to 
figure a) at different magnetic offset fields. 



eral, the effect of COM-REL motion coupling can have a 
significant impact on the stability of ID quantum gases 
and should therefore be considered in current ID exper- 
iments. 
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SUPPLEMENTAL MATERIAL 



Fitting the trapping potential 



Measurement of the single particle excitation 
frequencies in the trap 



We excite individual atoms in the ground state of the 
trap to higher trap states by periodically modulating the 
center position or the confinement length of the trap with 
a certain frequency uj. The transition frequencies are 
determined from fits to the excitation spectra (see [34]). 
The results are given in Table [il] 



h.o. number 


frequency 


(n x ,n y ,n z ) 


cj/2tt [kHz] 


(0,0,1) 


1.486(0.011) 


(0,0,2) 


2.985(0.010) 


(0,0,4) 


2.897(0.020) 


(1,0,0) 


13.96(0.08) 


(0,1,0) 


14.82(0.09) 


(2,0,0) 


26.43(0.27) 


(0,2,0) 


28.26(0.25) 



Table II. Transition frequencies in the trap of atoms excited 
in the trap. We denote the transitions by the corresponding 
quantum number of a harmonic oscillator. 



Determination of the elastic CIR 

To determine the elastic CIR we only consider one 
resonance following p~3j [14] although we have found an 
anisotropy in our system. We calculate the position of 
the elastic CIR from the transversal confinement length 



d i 



(8) 



To determine the mean trap frequency uj± we fit the 
Gaussian shape of the optical beam which creates the 
trapping potential to the measured transition frequencies 
(see [34 ). From the two profiles in radial direction we 
deduce the corresponding transition frequencies between 
the lowest single particle bound states in the gaussian 
potential and calculate the mean frequency 



2 (w ( i,o,o)fit + "(o,i,o)fit)=27r x (14.22 ± 0.35) kHz . 

(9) 

Using equation [8] as the input of the theory of [31] and 
the calibration a(B) of [30] we calculate the position of 
the elastic CIR. The propagated error of the mean trap 
frequency serves as the systematic uncertainty of the po- 
sition of the CIR. 



In an ideal and simplified case, the shape of the exter- 
nal potential is determined by the profile of a Gaussian 
beam which results in the potential 



V(r) = V 



w(z) 



exp 



-2(x 2 +y 2 ) 
w 2 (z) 



(10) 



where w(z) = w )J 1 + [j^j 



and wq = w(0) is the waist 
- is the Rayleigh range. 



size of the beam and = ™' 
The measurement of the excitation energies, Table [TTJ 
reveals an ellipticity of uJ X 2-o/^y2-o — 1-07 i n radial di- 
rection which cannot be justified by the beam potential 



Eq. (10). Hence, for the theoretical description, the po- 



tential Eq. (10) is expanded in a Taylor series around the 
origin up to second order which results in a harmonic 
potential of the form 



V(r) 



h 

2 



2 -2 
Uj3 ■ 



(11) 



To address the ellipticity in the radial direction, differ- 
ent harmonic ujj frequencies for the x and y direction are 
chosen. In the harmonic approximation, however, the 
coupling between COM and REL degrees of freedom van- 
ishes. This is surly not the case for the finite, anharmonic 
experimental potential. Hence, for the model potential 
higher order terms are introduced to include COM- REL 
coupling. It has been demonstrated that sextic poten- 
tials, i.e. expansions of a sin 2 optical lattice potential up 
to order six 



V(r) 



E 45 

=x,y,z 



-Vk 4 n 4 



(12) 



are well suited to describe single well potentials including 
COM-REL coupling [25]. 

The resulting fitting parameters (given in atomic units) 
are k x = 9.23 x 10~ 5 a.u., k v = 8.77 x 10~ 5 a.u., k z = 



1.58 x 10- 5 a.u., V x = 6.82 x 10~ 
a.u. and V~ = 2.33 x 10" 12 a.u. 



a.u., K = 6.58x10" 



COM-REL coupling involving the state |^ (b) ^ 2 ,o,o) 

The measurement of coherent coupling between the 
repulsive non-bound state and the COM-excited bound 
state \^ b ^ $2,0,0 ) nas been performed in a similar way as 
the one presented in the main text with the only differ- 
ence that an additional linear potential had been applied 
in z-direction (see Fig{4]). The additional potential reads 

Min = IhnB'z (13) 

where \i m is the Bohr magneton and B' = 18.92 G/cm 
a magnetic field gradient. The additional potential 
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slightly shifts the energy of the non- 
1-0(0) ^0,0,0) and thus also shifts the 
REL coupling resonances. Yet, the 
pling strength should be negligible, 
results for the Rabi frequency and 
in Fig. [4] for a comparison with the 
which has been performed without 
tial in z-direction. 
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the width presented 
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Figure 4. Dynamic of the COM-REL motion coupling, (a) 
Oscillation between the non-bound state and the COM ex- 
cited molecular state \ip^ $2,0,0)- (b) Maximum fraction of 
molecules extracted from measurements similar to a) at dif- 
ferent magnetic fields. 



